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ABSTRACT 



Context. The CoRoT satellite has made available high precision photometric observations of a large number of stars of different 
spectral types. Continuous photometric time series allow the characterization of stellar microvariability in a systematic way. 
Aims. We determine an index indicating the level of activity, derived from photometric data, for a large sample of stars with different 
color temperatures. We also assess to what extent this index can be related to an estimated Rossby number for stars whose rotation 
period can be estimated from the light curve. We also estimate a characteristic lifetime of the surface heterogeneities causing the 
variability of selected light curves. 

Methods. Most work on the photometric impact of stellar activity has been based either on measuring the variance of stellar 
microvariability, or on fitting the light curve. Following similar research for the Sun, our work is based on the Fourier analysis 
of stellar light curves. We have analyzed the Fourier power spectra of 430 selected light curves obtained by CoRoT during three 
observation runs. The low-frequency contribution of the stellar variability is modelled by a "generalized semi-lorentzian" profile. 
An activity index is derived from the fitted amplitude and width of the semi-lorentzian model. Some of the Fourier spectra exhibit 
a rotational modulation which enables the determination of the rotation period. In addition, a convective turnover time is derived 
from a grid of stellar models, so that a Rossby number can be estimated for this subsample of stars. A characteristic lifetime of the 
phenomena causing the observed power at low frequency is assessed from the fitted model of the power spectrum and is compared to 
the rotation period. 

Results. Higher values of the microvariability index are observed among the coolest stars from our sample. 28 light curves show a 
clear rotational modulation. The rotation periods derived from the observed low-frequency peaks in the Fourier spectrum decrease 
with the color temperature, in accordance with previous observations. The estimated Rossby number of most of the observed stars 
with a rotational modulation is less than 1, generally accepted as a critical value under which stars are expected to be active. The 
activity index, computed from the Fourier spectrum, decreases with increasing Rossby number. The quality of the CoRoT data 
enables the determination of the characteristic lifetime of active structures. It is shown to increase with the rotation period, but some 
scatter could arise from different surface heterogeneities (spot and faculae for example). 

Key words, stellar microvariability, stellar activity 



1. Introduction 



Since iKronl d 19471 [l950) suggested that spots on the surface of 
stars other than the Sun could explain their photometric vari- 
ability, much evidence for stellar activity has been collected. 
Activity refers here to a complex set of phenomena, including 
spots, faculae, eruptions, which have been observed on the Sun 
and are now suspected to exist on other stars. Activity is gener- 
ally considered to have a magnetic origin. However, the physical 
mechanisms that explain the observed active structures and their 
temporal evolution are still poorly understood. Observations of 
the Ca emission lines allowed the identification of stars with 
chromospheric activity, some of them having an activity cy- 
cle (IWilsonl lT978). Some time series of Ca lines fluxes show a 
rotational mod ulation, which is u s eful for determining the ro- 
tation period ( Vaug han et al.l [198 1 : Baliuna s et al.ll 19831) . as for 
the Sun. Co-rotating photospheric heterogeneities, such as stel- 
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lar spots, can also result in broad-band photometric microvari- 
ability. Large enough spots may be detected through the ro- 
tational modulation they induce in the star light curve. Due to 
their intrinsic temporal evolution, the surface heterogeneities 
are also expected to produce a power excess at low frequency 
in a Fourier sp ectrum of a velocity or photometric time series 
dHarvevl [ T9 85). These heterogeneities can be considered as a 
random process with some memory. They give rise to an expo- 
nentially decaying autocorrelation in time, yielding a Lorentzian 
shape in the power spectrum. In solar observations, at low fre- 
quency, several processes contribut e to the power spectr um: su- 
pergranulation (in velo city series , Jimenez et al. 1988), bright 
points (in photometry, Har vey etal.1 1993*1) or granulation (in 
both). As it is the process with the longest time scales, mag- 
netic activity seen through the signature of transiting spots on 
the solar surface in photometric series is detected at the low- 
est frequencies in the power spectrum. In solar-type stars, we 
expect the active structure s to be the phenomena with the 
longest time scales. Mosser et alj (12009) have shown that, in 
the case of several solar-type stars, the microvariability at 
the lowest frequencies is well modelled by the transit of stel- 
lar spots. In addition, simulation of the granulation and com- 
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parison with observations in the case of one of these solar- 
type stars (HD49933) sho ws that granulation has time scales 
shorter than that of spots (Ludwig et al. 2009), as for the Sun. 
Modeling of the spectrum with Harvey profiles is commonly 
used in the seismic analysis and shows a sim ilar shape to that 
of the Sun dMichel et alJl2009l:lBenomar et alj [2009). It is not 
excluded that other physical processes, such as large convec- 
tive cells for instance, which are suspected to exist in giant 
stars, could have the longest times. However, we selected our 
targets among stars that are identified as main-sequence. 
Stellar activity has been raised as an obstacle to th e detec - 



tion of exo-planets around cool stars ( e.g. Aigrain et al.l ((2004), 
lLagrange et all d2010h . iMeunier et all d2010l) ). The impact of 
stellar activity on exoplanet detectability has been estimated 
based on the solar case and on scaling laws (Ai grain et alj 
120041) . Simulations have given estimates of the "stellar ac- 
tivity noise" that strongly depend on spectral types. Actual 
light curves collected by the CoRoT (Convection, Rotation 
and planetary Transits) now provide us new insights on the 
microvariability that can be attribut ed to stellar activity. 

As stressed bv lRibarik et al. (2003), CoRoT gives the oppor- 
tunity to analyze a large amount of photometric data related to 
spotted stars. While ground-based observations have a limited 
precision and a strongly constrained observational window, the 
space-based CoRoT photometric instrument has collected nu- 
merous high-precision and continuous light curves. It has been 
designed to study several bright seismic targets and to search 
for exoplanets around a large amount of fainter stars. Seismic 
and exoplanet targets belong to different spectral types and lu- 
minosity classes. Exoplanet targets have however been mainly 
selected from F to M spectral types. The observed noise from 
an exoplanet search point of view has been measured and has 
been show n to be from two to three times the expected pho- 
ton noise (Aigrain et all 120091) . We have performed a more 
detailed analysis of the stellar microvariability based on the 
Fourier power spectrum, which allows us to discriminate 
several phenomena with different time scales. The purpose of 
this work is to show, based on a still limited but homogenous 
sample, that the available light curves can deepen our knowl- 
edge of stellar activity as the time scales of magnetic structures 
provide hints to understand the underlying mechanisms. 

There are at least two approaches to process the thousands 
of available light curves. The first is based on direc tly modeling 
a t ime series, as sugg ested by Ribarik et al] (l2003|) and applied 
bvlLanza et all d2009l) . iMosser et all d2009h and lFrohlich et alj 
J20O9ir Models of light curves are based on various assumptions 
regarding the number of spots, their physical properties, differ- 
ential rotation and limb darkening. There is no unique solution, 
but some "reasonable" hypotheses can give interesting results 
for a case by case analysis. A second a pproach is bas ed on the 
Fourier power spectrum of a time series (Harvev 1985). It can be 
more reliable in identifying different features in a given signal 
as it can discriminate different time scales and quantify the 
lowest frequency contributions to the stellar microvariabil- 
ity. In addition, it allows the identification, in some cases, of a 
rotational modulation. This second approach is adopted in the 
present work. Our work is based on a 430-star sample, which is 
large enough to obtain several interesting results regarding stel- 
lar activity. Here we develop and test an analysis method on the 
present sample in orderto determine physical constraints on stel- 
lar activity. 

In Section |2 we present our initial sample of 430 stars and 
describe our analysis method of the Fourier power spectra of 
the selected light curves. In Section|3] we first summarize some 



useful results related to the solar microvariability. We then define 
our stellar activity index computed from the Fourier spec- 
trum. The distribution of our activity index against the color 
temperature is analyzed. In Section |4] we present our findings 
related to 28 light curves, out of the 430-star sample, show- 
ing a rotational modulation. In Section|5] the Rossby number for 
stars showing a rotational modulation is derived. The correlation 
between the observed activity index and the estimated Rossby 
number is investigated. In Section [6] we analyze our results re- 
garding the characteristic evolution times of microvariability. In 
Section|7] we conclude and discuss the opportunity to go further 
into a wider analysis of the CoRoT data. 



2. The sample of CoRoT light curves 

2. 1 . Selection criteria 

Our sampl43 was initially made of the first 30 seismic targets 
observed by CoRoT during three observation "runs" (IRaOl - 
January/April 2007- SRaOl -April/May 2007- and LRcOl - 
May/October 2007) and of the 400 brightest exoplanet targets 
observed during the initial run (IRaOl). The apparent V magni- 
tudes of the seismic targets range from 5.45 to 9.48, and between 
12 and 13 for the selected exoplanet targets. Our sample includes 
stars with spectral types from B to M, but most of the exoplanet 
targets are redder than A. Most of the selected light curves are 
60-days long. However, some of the seismic targets light curves 
are limited to 30 days and a few are 150 days long. Seismic 
targets are observed in a broad-band white-light flux. For exo- 
planet targets, the flux is the sum of three broad-band chromatic 
fluxes, equivalent to a white-light flux. There is no bolometric 
correction. Broadband photometric data are collected with two 
sampling times (32 or 512 s) depending on the kind of target. 
The available photometric time series have a very high preci- 
sion that allows, after correcting for an instrumental linear trend 
with time, a systematic study of stellar microv ariability. More 
details on the CoRoT time series are given by Auvergne et al.l 
(120091) . 

2.2. Available stellar information 

Apparent B, V, R and I magnitudes were determined during 
ground-based p hotometric c ampaigns in preparation for the 
CoRoT mission dDeleuil et al.l 120091) . A color temperature was 
derived from those magnitudes, which w ere not corrected for in- 
terstellar extinction dDeleuil et al.ll2009l) . We checked that there 
is a monotonic relationship between this color temperature and 
the observed B-V index of the selected stars. In the next sections, 
we use this temperature as a proxy for the stellar spectral type. 

2.3. Analysis procedure 

Each light curve has been analyzed through a 4-step process. 
The first step starts with a visual inspection of the light curve. It 
is performed in order to identify any significant discontinuities 
of the recorded flux which could be caused by a high-energy 
particle on the photometric detector. As the resulting disconti- 
nuities do not show any systematic time scale, no automatic 
procedure worked well enough to remove their impact in the 
time series or in the Fourier power spectrum. Consequently, 
a visual inspection has been performed (for more details 



1 All data are available from the public CoRot archive at: 
idoc-corot . ias . u-psud . fr 



J.C. Hulot et al.: Stellar activity from CoRoT data 



3 



C O/B A F C- K/M 

Spectral type 



Fig. 1. The distribution of the selected stars by spectral type. 



on these discontinuities, see Auverg ne et al J 12009). The light 
curves showing such discontinuities followed by a slow recov- 
ery are excluded from the next steps of our analysis. Exceptions 
are accepted if the discontinuity occurs close to the extremities 
of the time-series, so that the truncated time-series, excluding 
the discontinuity, is long enough to be included in the further 
analysis. All of the seismic targets went through the first step, 
but 67 exoplanet targets had to be excluded. A linear variation 
was removed with an ordinary least square linear fit. The ob- 
served temporal drift is mainly attributed to instrumental aging. 
We checked that a model with a second or third order polynomial 
drift does not significantly change the Fourier power spectrum. 
Missing observations are replaced by linearly interpolated val- 
ues. They account for close to 10 % and mainly result from the 
crossing of the South-Atlantic Anomaly by the CoRoT space- 
craft. Their signature appears in the Fourier spectrum at the 
orbital frequency and harmonics. The seismic targets include 
4 stars that were found out to be binaries. They were excluded 
from our analysis. At the end of this first step, our sample in- 
cludes 359 stars (Table [TJ which are mainly distributed within 
the A to G spectral types (Figure[TJ- 

In a second step, we derive the Fourier power spectrum from 
a uniformly sampled time series. The third step consists in com- 
puting several indexes to characterize the low-frequency power 
spectrum. Narrow spectral peaks are searched for and their fre- 
quencies are compared in order to identify the rotation period 
(see Section 3J. The fourth step is designed to fit a general- 
ized semi-lorentzian model to the low frequency spectrum. The 
power spectra of 58 stars, among the hottest ones from our sam- 
ple, were however too flat to be fitted to our model. The corre- 
sponding stars were probably not active or had too low activity 
to be detected. The number of stars surviving the different steps 
of our analysis is shown in Table Q] 

3. Low-frequency microvariability 

3.1. The solar case 

The solar microvariability has been extensi vely analyzed based 
on the SOHO/VIRGO photometric data (lAigrain et alj I2004 
ISeleznvov et alJ [2003) before being extrapolated to other 
stars (using scaling laws to describe the rotation, the ampli- 
tude and characteristic evolution time of activity, in order to 
simul ate an "activity noise " impact on exo-planet detectabil- 
ity by Aigrain et al.l d2004l) . VIRGO data can be compared to 
CoRoT data, as they are also intensity measurements. However, 
the effects of the instruments (bandwith) are different and should 
be taken into account. VIRGO data are based on 3 narrow band- 
widths, whereas Co RoT relies on one wide band. As shown by 
Aigr ain et al.ld2004 |). VIRGO power spectra show several contri- 



butions which are attributed to solar activity, super-granulation 
and granulation. The first contribution is observed at the lowest 
frequencies and shows a clear variability throu gh the solar cy- 
cle. A lorentzian profile was first suggested by iHarvevl (fl985) 
to model the low-frequency background in such a power spec- 
trum. Each contribution to the power spectrum is actually better 
described by a 'generalized' lorentzian profile, the total profile 
being the sum of each contribution: 



p(f) = 2 



- Htt)" 



+ B 



(1) 



P(f) is the power density in ppm^Hz 1 , with ppm for part- 
per-million. a, equals 2 when the phenomenon causing the flux 
variation exponentially decays, but different values can result 
from more complex temporal evolution. A, is the amplitude of 
the i-th contribution and 1/Af, is its characteristic time. B corre- 
sponds to an additional white noise. It is important to stress that, 
in the solar case, t he low frequency co ntribution has a 10-day 
characteristic time (Aigrain et al. 2004). It is then significantly 
shorter than the evolution time of a solar active region, generally 
estimated around 2 months, and should be interpreted as a char- 
acteristic evolution time of individual spots and faculae forming 
an active region. 

3.2. A simplified model of the low-frequency stellar 
microvariability 

A variance derived directly from a light curve could capture 
different superposed phenomena. On the contrary, a Fourier 
analysis allows the identification and quantification of the low- 
frequency background contribution. What is observed in the so- 
lar case is expected to be seen in some stellar light curves if 
the photometric precision is high enough. The low-frequency 
part (below 100 //Hz) of the power spectrum may result from 
several contributions that must be disentangled. A first contri- 
bution to the power excess at low frequency results from in- 
strumental effects. All our light curves are corrected for a lin- 
ear trend and the CoRoT stability is confirmed by a significant 
number of flat power spectra, mainly observed for A stars. A 
second contribution is from the so-called "hot pixels" that re- 
sult from incident high-energy particles on the photometric de- 
tector. We recall that the light curves concerned are excluded 
from our analysis through a visual inspection, so that this contri- 
bution should not be observed in our sample. A third contribu- 
tion to the power spectrum relates to intense and narrow peaks 
that are due either to long-period oscillations or to a rotational 
modulation. The fourth contribution is the one that we try to 
identify and to model. It corresponds to a low-frequency back- 
ground and is well described by a continuous profile centered on 
zero frequency (see Eq.Q}. It varies significantly from one light 
curve to another. Our sample shows that an A-star spectrum is 
more often flat, while a G-star spectrum shows a higher power 
at low frequencies. We suppose that late-type main sequence 
stars show surface heterogeneities with time scales that fol- 
low the same hierarchy as in the solar case, that is activity 
has a longer time scale compared to supergranulation and 
granulation. This relies on previo us observations: firs t, as al- 
ready mentioned, the results of Mosser et al. (2009) about 
spot modeling; second, in the case of HD49933, the power 
spectrum expected from a granulation model has been com- 
pared to the actual power spectrum observed with Co RoT 
and show the expected hierarchy of time scales (Ludwi g^t alJ 
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step 


number of 


number of 


total number 




seismic targets 


exoplanet targets 


of stars 


initial selection 


30 


400 


430 


after excluding 'hot pixels' and binaries 


26 


333 


359 


fitted Fourier power spectrum 


25 


276 


301 


with rotational modulation 


6 


22 


28 


with estimated Ro 


6 


21 


27 



Table 1. Number of selected light curves through the different steps 



l2009h . However, we cannot exclude that a giant star is in- 
cluded in our sample and that its low-frequency microvari- 
ability comes from another phenomenon than magnetic ac- 
tivity. Such a case is likely to happen within some red giants 
whose surface could show large convective cells with long 
evolution time scales. Our sample of stars has been selected 
from main-sequence stars from the CoRoT database but mis- 
classifications cannot be excluded. 

Considering the power spectrum at frequencies lower than 
lOOyuHz, we exclude the peaks that may result from a rotational 
modulation or from long-period oscillations. This 'cleaned' 
power spectrum is then fitted to a single generalized semi- 
lorentzian distribution to which we add a constant free param- 
eter: 

P(f) = — T7^ +S (2) 

A is the amplitude of the lower-frequency contribution to the 
spectrum and Af is related to its characterist ic time. It cap- 
tures as shown in the solar case (lHarvevl [l985). the stellar ac- 
tivity. B captures all the other contributions, including granula- 
tion, and not only the additional white noise as in Eq.[TJ 

3.3. The low-frequency background index 

We derive a microvariability index from the generalized semi- 
lorentzian model (Eq. We call it a "low-frequency back- 
ground index" (LFBI). It is defined as Iif = AAf and it is ex- 
pressed in ppm 2 . It represents the integral of the power density 
of the model which describes the microvariability disentangled 
from other contributions. A possible bias in the computation 
of the index corresponds to the case of a slowly rotating star 
whose light curve has an unresolved rotational modulation. 
In such a case, the rotational modulation peak cannot be re- 
moved from the power spectrum: it affects the very first bins 
of the spectrum and may lead to an over-estimated LFBI. 
In order to use the solar case as a reference when consider- 
ing CoRoT light curves, we have analyzed some VIRGO data. 
These data need to be corr ected for the instru ment bandwith 
(see Section [3~TT i. following iMichel eta l. (2009). The correction 
factor to VIRGO data to be compared with CoRoT data is 
estimated to be 1 .4. The results obtained from the VIRGO solar 
data at different times of the solar cycle indicates that AAf is 
a representative index for activity. In the solar case, while Iif is 
close to 1 .4 10 5 ppm 2 at the maximum of activity, it is nearly 100 
times lower at its minimum. Correspondingly, A varies from 
1.5 10 5 ppm 2 ^Hz 1 at maximum to nearly 100 times lower at 
minimum, while Af remains close to 1, which corresponds to 
nearly 10 days. 

In order to check the robustness of our index Iif, two addi- 
tional indices were computed, consisting of the direct integra- 
tion of the Fourier power spectrum, respectively from up to 
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Fig. 2. HD 49933 light curve. The observed signal has been cor- 
rected for a linear trend and converted into a relative flux varia- 
tion and presented in ppm. 



10 ;uHz and from up to 100 yuHz. They were found to be highly 
correlated with Iff, but we recall that these latter indices take 
into account all of the power at low frequency, whereas Ijji does 
not take into account the power due to pulsation peaks, peaks 
due to rotational modulation or the constant term capturing other 
phenomena and noise. 

3.4. A study case: HD 49933 

HD 49933 is a bright solar-like star for which several parame- 
ters are now well known from ground-based observations. It is 
an F5 star with a visual magnitude 5. 77. Its effective temper- 
ature is between 6500 K and 67 80 K i lRvabchikova et alj [2009: 
iBruntt et all [20081: iBrunttl 12009b . HD 49933 was observed by 
CoRoT during its 60-day initial observation run with a 32s sam- 
pling time. The de trended stellar intensity shows some signif- 
icant variations with a ~1000 ppm amplitude (see Fig. 0. We 
note that the VIRGO photometric data show a peak-to-peak vari- 
ation which is also close to 1000 ppm at the maximum of solar 
activity. The power spectrum of HD 49933 light curve confirms 
a power excess at low frequency (see Fig. [3). 

Following the method presented above, we obtained the 
following fitted parameters (see Fig. |4j, using a MLEfit- 
ting algorithm described by Anoourchaux et al. (2008): 
A = 1.010 4 ppmVHz- 1 , B=1.7ppm 2 J uHz- 1 , A/ = 2/zHz and 
or = 2.4. A is about 15 times lower than the observed value for 
the Sun at its activity maximum, but less than 10 times larger 
than the same parameter at the Sun's minimum of activity. The 
HD 49933 activity index is 2. 10 4 ppm 2 . Observations on a long 
time scale wo uld be necessary to get a better knowledge of its 
activity cycle (iGarcia et al ] l2010HWilsoiir9 78). 
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Fig. 3. Power spectrum of HD 49933 computed from a 60-day 
time series after it was corrected for a linear trend which cap- 
tures some instrumental effects. The power is presented here on 
a logarithmic scale that allows us to see the rotational modu- 
lation peak at 3.5 + 0.2 fiHz and its third harmonics. . 



Fig. 5. The fraction of most active stars, defined as stars with a 
LFBI higher than a given threshold (see text). Here, the thresh- 
old was defined so that the most active stars account for 20 % 
of the total sample. The fraction of stars with an index higher 
than the threshold is larger for the coolest stars. This general 
trend is not sensitive to the threshold value. 




Frequency (/iHz) 



Fig. 4. A 'generalized' lorentzian model (with a uniform back- 
ground) is fitted (solid line) to the power spectrum of HD 49933 
(crosses), from which the rotational modulation peaks have been 
removed. 



3.5. Color temperature and microvariability 

The LFBI was computed for the 301 stars of the sample for 
which a fit could be performed. As was noted in section 2.3, 
58 power spectra could not be correctly fit. All of them were 
"flat" spectra with a low power density at low frequency. 
Such spectra confirm the instrument stability. The corre- 
sponding stars are among the hottest and do not show a 
photometric signature of surface heterogeneities. The activ- 
ity indexes do not show a clear dependence on color tempera- 
ture. However, our sample is not uniformly distributed in terms 
of temperature: F stars are over-represented. In order to avoid 
this bias, we sort our sample by increasing temperatures, and 
divide it into nine sub-samples, each of them including nearly 
the same number of stars. We determine a threshold value of the 
LFBI defined so that 20% of the full sample has an index higher 
than the threshold. We then determine, within each temperature 
sub-sample, the fraction of stars with a LFBI higher than the 
threshold and labelled by the fraction of most active stars. The 
fraction of most active stars appears to be clearly related to the 



temperature, being higher in the low temperature sub-samples 
(see Fig. [5]l. As the 20% threshold is arbitrary, we checked that 
the same result is observed when the threshold is determined so 
that we select respectively the 10%, 30%, 40% or 50% most ac- 
tive stars from the whole sample. The higher fraction of active 
stars among the G to K spectral typ es is consistent with previ- 
ous observations (Berdyugina 2005). According to our results, 
cool stars are more likely to show an enhanced microvariability. 



4. Rotational signature 

4.1. Rotational narrow peaks 

Some of the analyzed power spectra show super-imposed narrow 
peaks in addition to the low frequency background. Such peaks 
are attributed to the transit of surface heterogeneities which 
are large and contrasted enough to cause a periodic intensity 
modulation over a significant fraction of the observation (see 
Moss er et al.l l2009). A Fourier analysis shows that a single co- 
rotating permanent spot results in a strong peak and weaker even 
harmonics in the Fourier power spectrum of its photometric flux 
(Clarke 2003). Indeed, we often observed harmonics above the 
peak corresponding to the rotation period. We checked with sim- 
ulated light curves that a more complicated distribution of spots 
or, more generally of active areas, gives the expected peaks, with 
even and odd harmonics. These peaks are thus considered as the 
signature of activity and rotation. Because the observation time 
is limited, the frequency resolution might be too low to detect 
such peaks in slowly-rotating stars. As previously noted, the 
LFBI of such stars could be over-estimated. 



4.2. An example of rotational modulation : HD 49933 

In order to identify a rotational modulation in the CoRoT light 
curves, we first focused on the brightest stars which were se- 
lected as seismological targets, and chose HD 49933 among 
them. While p-modes were detected in this star, a rotational 
modulation of the observed spect ral lines was attributed to one 
or more spots (Moss er et al . 2005). In order to illustrate how we 
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spectral 


number 


out of which with 


proportion with 


type 


of stars 


detected rotational 


detected rotational 




step 1 


modulation 


modulation 


O/B 


9 


o 


0% 


A 


129 


2 


2% 


F 


140 


13 


9% 


G 


57 


12 


21 % 


K/M 


24 


1 


4% 


Total 


359 


28 


8 % 



Table 2. Number of selected light curves with a detected rota- 
tional modulation signature. The spectral type is derived from 
a color temperature determined from ground-based photometric 
observations. 



carried out our analysis on the full sample of stars, we present in 
more detail our results for this star. 

The power spectrum in the range [0,20] fjHz shows a strong 
peak and weaker ones (see Fig. [3}. The main peak is centered at 
3.5 + 0.2 yuHz which corresponds to a rotation period of 3.4 + 0.2 
days and is co nfirmed by the observed p modes oscillations. 
dAppourchaux et aLll2008l) . A second peak is observed at twice 
this frequency and a third one at a frequency four times larger. 
The second harmonic is ten times weaker than the fundamental 
but twice as strong as the fourth harmonic. We see neither a 
third harmonic, nor a fifth. Such a result would be obtained with 
a single and long-lived co-rotating active region (Clarke 2003) 
and is also in agree ment with the analysis of the light curve by 
Moss er et al ] (120091) . 

4.3. Rotational modulation within our sample 

There may be active stars whose power spectra do not show ro- 
tational peaks. A first reason would be an axis of rotation too 
close to the line of sight. A second reason is a rotation period 
that is too long to be identified in a too low-frequency reso- 
lution Fourier power spectrum. A third reason could be that 
stars having spots evolving on time scales shorter than the 
rotation period do not show a clear rotational modulation. 
Indeed, our analysis in terms of LFBI shows cool stars with a 
high excess power at low frequency, but without any clear ro- 
tational modulation. However, 6 stars (including HD 49933) out 
of 30 from the seismological sub-sample and 22 stars from the 
exo-planets sub-sample show a rotational modulation. The frac- 
tions of identified rotational modulations in the two sub-samples 
may be biased by their magnitude distributions. The asteroseis- 
mology targets are significantly brighter than the exoplanet ones, 
so that we detect weakly active stars within the former, but only 
strongly active ones within the latter. We identified two G stars 
with a period too long (respectively more than 23 days and more 
than 29 days) to be confirmed with a 60-day time series. We de- 
cided not to include these two stars in our sub-sample of stars 
with a rotational modulation. 

The rotational modulation appears to be significantly more 
frequent among G and, to a lesser extent, F spectral types (see 
Table|2|. The low proportion of K and M stars with an observed 
rotational modulation seems to be at varianc e with the expecte d 
activity of K and M stars (see for instance Berdvugina 2 0051) . 
However, our sample of K and M stars is too limited and their 
rotation periods are likely to be too long for a rotational modula- 
tion to be observed in a 6 0-day photometric time-series. Indeed, 
Kiraga & Stepien (2007) identify several M stars with 30-day 
and longer rotation periods. 
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Fig. 6. Rotation periods versus the color temperature of the 28 
stars whose light curves show a rotational modulation. 



The rotation periods of stars with an activity signature are 
derived from the rotational peaks. In most of the observed cases, 
the frequency resolution is too low to determine the width of the 
peak that is attributed to the rotational modulation. In two cases, 
we can only determine a lower estimate of the rotation period, 
so that the corresponding stars are excluded from the next steps 
of our analysis. When plotted against the color temperature (see 
Fig- IS), the rotation period shows a clearly decreasing trend ver- 
sus the color te mperature . This result is consisten t with previous 
observations dGravlll982t lLockwood etaDT l984). A significant 
scatter results from additional parameters which are thought to 
impact the rotation period (the age for example). 



5. Rossby number 

5.1. A determinant of stellar activity? 

To exhibit solar-like activity, a star needs to have a convective 
envelope. According to our present understanding of stel- 
lar interiors, stars redder (cooler) than F2 have a convective 
envelope whose thickness decreases withan increasing sur- 
face temperature. A dynamo effect is expected to take place 
when the dynamo number, which compares advective to dif- 
fusive effects on the magnetic fields, is larger than 1. Under 
some very simplified assumptions, a dynamo number can be 
expressed as a decreasing function of the Rossby number 
(lParkerlll979b . 

The Rossby number is the ratio between the surface ro- 
tation period and the convectiv e turnover tim e at the bot- 
tom of the convective envelope. Gilliland (1985) actually ob- 
served that a chromospheric activity index, based on the 
CallH and K emission fluxes, is correlated with an esti- 
mated Rossby number within a 41-star sample. The convec- 
tive turnover time is expected to be a decreasing function of 
the effective temperature. A rough relationship between the 
B-V index and the rotation period can be inferred for main- 
sequence stars from a statistical study of a large population 
of stars (Gravl il982h . On average, a cooler star rotates more 
slowly than a hotter one. However, the rotation period and 
the convective turnover time are likely to show significant 
scatter within a given spectral type, so that the B-V index is 
not enough to determine if a star is active or not. The Rossby 
number is expected to be relevant as an activity index. 
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From a 277-star sample for which groun d-based photo- 
metric observations are available. lHall £1991) concluded that 
the Ros s by nu mber of active stars must be lower than 2/3. 
Stepien (1994) showed that an activity index is well corre- 
lated with the Rossby number for a limited color class of 
main-sequence stars, but he also stressed that, for the other 
stars, this dimensionless number does not explain activity 
better than the rotation period. As noted bv lHalll ( 119911) . the 
conditions for a star to be active have generally been ad- 
dressed focusing on samples in which (very) active stars are 
over-weighted. A more systematic study of a larger sample 
of stars, as made possible with CoRoT, spread over the HR 
diagram could avoid such bias. 



5.2. Estimating a Rossby number 

The Rossby number is defined as the ratio of the rotation pe- 
riod to the convective turnover time at the bottom of the convec- 
tive envelope (jcv)- The rotation period is derived from the fre- 
quency of t he main peak in th e Four ier power spectrum. Looking 
at M stars, Kiraga & Stepien] (120071) suggested an empirical de- 
termination of the convective turnover time based on a linear 
relation between the logarithm of an X ray-fl ux index, not avail- 
able for our targets, and the rotation period. iNoves et al. ( 1984) 
used a scaling relation between a color index and a convective 
turnover time. 

Here, we prefer to estimate this turnover time using a grid 

of stellar models with 1 to 2 solar masses, for a star whose 

effective temperature is presumed to be known + 100 K. The 

global ID models are obtained with the CESAM code (ver- 

i ii 1 

sion 4, see Morel & Lebreton 2008) assuming standard physics. 

Convection is described according to iBohm- Vitensel (U 95 8 ) ' s lo- 
cal mixing-length theory of convection (MLT) with a mixing- 
length A = a c H p , where H p is the pressure scale height 
and a c is the mixing-length parameter. The calibration of the 
associated solar model gives a c = 1.62. Turbulent pressure 
and microscopic diffusion are not included. All models have 
a solar iron-to-hydrogen abundan ce and the chemica l mixture 
of the heavy ele ments of iGrevesse & Noelsl d 1993b . We use 
OPAL opacitie s (Iglesias & Rogers! Il996l) extended with the 
I Alexander & Ferguson! (11994 data for T < 10 4 K , both sets 
of data being given for the Grevesse & Noels (1993) solar mix- 
ture. Finally, the CEFF ( Christensen-Dalsgaard & Dappenl l992) 
equation of state is assumed. The convective time Tcv is then 
directly deduced from the mixing length A and the convective 
velocity u at one height scale from the base of the convective 
zone. There is indeed a difficulty in estimating a turnover time at 
the bottom of the convective envelope where the velocity is zero, 
which is why we estimate it at a height H p above the bottom. 
The Rossby number and the low-frequency background index 
of the Sun are useful references. The solar convective turnover 
time is computed with the same grid of stellar interior model . We 
obtain a Rossby number close to 1 . The LFBI for the Sun is de- 
rived from VIRGO data for a 180-day active period. We obtain 
log(LFBI) ~ 5, which places it among very moderately active 
stars. 



5.3. Rossby number and microvariability index 

In our sample, one star was expected to be fully convective so 
that our estimation concerns 27 stars from our 28 stars with a 
rotational modulation. Fig. [7] shows a clear anti-correlation be- 
tween the activity index and the Rossby number: the most active 



3' 



Rossby number 



Fig. 7. The low-frequency background, or microvariability, in- 
dex versus the Rossby number. The active Sun, i.e. at its maxi- 
mum, corresponds to the square, while the quiet Sun would have 
an index one hundred times lower (see section 3.3). 



stars having the lower estimated Rossby number. The observed 
scatter could result from various parameters, one of which 
is the stars metallicity. A similar relation was searched for 
between the rotation period and the activity index. No trend 
was observed. The Rossby number seems then to be a more 
useful predictor of the low-frequency microvariability than 
the rotation period. The Rossby numbers are limited by the 
observed rotational periods and by the color temperature of the 
selected stars. Based on the CoRoT photometric data from the 
initial run, our analysis method allows the determination of the 
rotation period over an interval from a few to nearly 30 days. 
The upper limit is quite restrictive for the coolest stars. It will be 
enhanced by using the CoRoT data from longer observation pe- 
riods. On the other hand, the spectral range studied corresponds 
to a wide range of convective turnover times (from 0.05 day to 
170 days). 



6. Characteristic evolution time of the surface 
heterogeneities 

6. 1 . Low frequency power density and characteristic 
evolution time 

While stellar activity has been mainly addressed in terms of ac- 
tivity index (in Ca II H and K flux, X-ray flux or broadband pho- 
tometric flux) and activity cycle, based on long observations, 
there is much less available data about the characteristic evolu- 
tion time of the surface heterogeneities. The parameter t — 1 / A/ 
determined when fitting the Fourier power spectrum to a gen- 
eralized semi-Lorentzian has been previously used to compute 
an activity index. The parameter itself gives interesting informa- 
tion about the characteristic evolution time of the surface hetero- 
geneities. 

When there are two different phenomena, for instance spots 
and faculae, or when spots evolve on a time scale that is shorter 
than the evolution time of active areas, the lorentzian contribu- 
tion may be significantly wider than the peaks resulting from the 
rotational modulation. Such a result is observed in most of the 
identified stars with rotational modulation. For instance, in the 
case of HD 49933, Af = 2juHz corresponds to a 5.5-day charac- 
teristic time. This value is slightly less than two rotation periods, 
whereas the width of the rotational modulation peaks is com- 
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Fig. 8. The characteristic time t against the rotation period, 
showing a clear increasing trend of t with P rot . The slow ro- 
tator with a rapid evolution time, indicated by a square, is a 
Ml star. 



parable to the frequency resolution of a 60-day time series. The 
rotational modulation could be caused by a phenomenon whose 
evolution time is quite long compared to the rotation period, for 
example faculae. 

6.2. The observed characteristic evolution times 

The characteristic time (r) is plotted against the rotation period 
(P) in Fig. [8] showing a clear increase of r with P. Several stars 
show a rotational modulation caused by a phenomenon whose 
characteristic time is very large compared to the rotation pe- 
riod. In our sample, very few stars show a rotational modula- 
tion caused by phenomena with rapid evolution, on a time-scale 
shorter than the rotation period. However, there is an interest- 
ing exception, a Ml star with a 11.6-day period and a 2.5-day 
characteristic evolution time. M stars have a deep convective en- 
velope or are fully convective. Even if they are generally slow 
rotators, they are expect ed to be very active, as confirmed by X- 
ray flux observations (Kiraga & Stepien 2007). A more exten- 
sive sample might allow us to show several different relation- 
ships between rotation period and characteristic evolution 
time. It might be explained by dilferent physical mechanisms 
linked to surface heterogeneities. 



7. Conclusion 

7.1. Summary and discussion of our results 

Our first goal was to test an analysis method designed to iden- 
tify and quantify an activity signature in the photometric data 
collected by CoRoT. Based on a homogeneous sample of more 
than 350 stars, we obtain original results regarding a contribu- 
tion to stellar microvariability from stellar activity. CoRoT, 
thanks to its sensitivity and stability, allows the investigation 
of weaker activity stars than from ground-based observa- 
tions. 

Our work shows that CoRoT data allow the detection and 
measurement in many cases of a low-frequency power ex- 
cess, which is interpreted as a photometric signature of stellar ac- 
tivity. After removing several contributions to the low-frequency 
power excess, we are able to fit the low-frequency background 
to a generalized semi-lorentzian model. A low-frequency back- 



ground index (LFBI) is then determined in a systematic way. 
Such an index is more reliable than a microvariability index di- 
rectly derived from the light curve, because it captures only the 
low-frequency contribution which is likely to result from stellar 
activity and excludes other possible sources of variability, such 
as low frequency pulsations. 

Despite some scatter, our activity index is correlated to the 
color temperature, the highest values being obtained for the 
coolest stars. This r esult is consistent with previous observa- 
tions, as reviewed by Berdvugina (2005) for example. Our sam- 
ple allows us to quantify more precisely this relation between the 
color temperature (or any parameter locating a star on the main- 
sequence) and its level of activity. However, a single-parameter 
analysis is likely to show significant scatter. Several factors, such 
as rotation or age and possibly metallicity through its influence 
on convection, may induce a scatter of an activity index against 
the temperature. Another possible cause for the observed scat- 
ter could be the cyclicity of stellar activity: the duration of 
our observations is likely to be shorter than the activity cy- 
cle. 

Our method allows us to quantify an activity index even if 
there is no rotational modulation. However a detected rotational 
modulation permits us to go further by estimating the Rossby 
number. We have seen a clear rotational modulation in the light 
curves of 28 out of the 430 stars initially included in our sample. 
There is a higher proportion of such stars within the G, and to 
a less extent, F spectral types. It could be explained by two bi- 
ases: hotter stars are less likely to be active and the cooler ones 
are expected to be slow rotators, so that their rotational modula- 
tion is not detected in a 60-day light curve. The most active stars 
have a low Rossby number, sometimes greater than unity but 
still anti-correlated with the activity index of the sta r. This re- 
sult is consistent with previous work (for instance lHalll ll 99 lh 
but our sample allows the exploration of larger values of the 
Rossby number. Some scatter is also observed, again as expected 
from a single-parameter analysis. The Rossby number is only a 
proxy for a dynamo number. The scaling relation between the 
two dimensionless numbers is based on the strong assumption 
regarding the differential rotation in the convective envelope that 
may be too approximate for some spectral classes. Moreover, 
the scatter of the activity index at a given Rossby number can 
be explained by adverse observing conditions (low activity dur- 
ing a cycle or axis of rotation close to the line of sight). In ad- 
dition, no clear relation is found between activity and rota- 
tion period in our sample, confirming the relevance of the 
Rossby number. However, we also investigated the variation 
of the characteristic evolution time, which is clearly corre- 
lated with the rotation period. 

Identifying active stars and sorting them according to their 
physical conditions is an interesting goal in itself. 

7.2. Perspectives 

Our sample covers a large part of the HR diagram, but is not un- 
biased. The exoplanet targets have been mainly selected among 
the cooler main-sequence stars. The available data however 
cover A to M spectral types. In any case, the study of a larger 
sample must take into account the non-uniform distribution in 
terms of color temperature. 

In order to study the correlation between an activity index 
and an estimated Rossby number, we need to determine the rota- 
tion period. Some power spectra do not show any low-frequency 
peaks, but have a high-amplitude low-frequency power excess. 
The corresponding stars are likely to be active. As was ex- 
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plained for M stars, a power excess at very low frequency can 
result from an unresolved modulation peak which corresponds 
to a slow rotation. When an activity index is plotted against 
the color temperature, the non-rotationally modulated but active 
stars are taken into account. Unresolved modulation peaks may 
over-estimate the corresponding activity index. To study the K 
and M stars more precisely, longer time-series are needed. More 
generally, a larger sample of stars is needed to describe more 
precisely the relationship between the characteristic evolution 
time, the rotation period and the color temperature (or any p a- 
rameter that locates a main-sequence star in the HR diagram). 
M stars are again particularly interesting as they are supposed to 
be fully convective and thus cannot be described with the same 
definition of the Rossby number. In slightly hotter stars, the con- 
vective turnover time is strongly sensitive to the color tempera- 
ture, since we consider stars with a shallow convective envelope, 
which are expected to show stee p variations of the convective 
turnover (Kiraga & Stepien 2007). 

CoRoT has now collected numerous longer time series and 
Kepler is collecting even longer time-series, both missions are 
thus very promising for the understanding of stellar magnetic 
activity. 
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